Reduced fidelity in topological quantum phase transitions 
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We study the reduced fidelity between local states of lattice systems exhibiting topological order. 
By exploiting mappings to spin models with classical order, we are able to analytically extract the 
scaling behavior of the reduced fidelity at the corresponding quantum phase transitions out of the 
topologically ordered phases. Our results suggest that the reduced fidelity, albeit being a local 
measure, generically serves as an accurate marker of a topological quantum phase transition. 
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Introduction — Electron correlations in condensed mat- 
ter systems sometimes produce topologically ordered 
phases where effects from local perturbations are ex- 
ponentially suppressed The most prominent exam- 
ples are the fifty or so observed fractional quantum Hall 
phases, with the topological order manifested in gapless 
edge states and excitations with fractional statistics. The 
fact that the ground state degeneracy in phases with non- 
Abelian statistics cannot be lifted by local perturbations 
lies at the heart of current proposals for topological quan- 
tum computation 

The insensitivity to local perturbations invalidates the 
use of a local order parameter to identify a quantum 
phase transition out of a topologically ordered phase. At- 
tempts to build a theory of topological quantum phase 
transitions (TQPTs) — replacing the Ginzburg-Landau 
symmetry-breaking paradigm — have instead borrowed 
concepts from quantum information theory, inparticular 
those of entanglement entropy Q and fidelity [ij , none of 
which require the construction of an order parameter. 

Fidelity measures the similarity between two quantum 
states, and, for pure states, is defined as the modulus of 
their overlap. Since the ground state changes rapidly at 
a quantum phase transition, one expects that the fidelity 
between two ground states that differ by a small change 
in the driving parameter should exhibit a sharp drop. 
This expectation has been confirmed in a number of case 
studies]!, including several TQPTs [l,i,[3. 

Suppose that one replaces the two ground states in a 
fidelity analysis by two states that also differ slightly in 
the driving parameter, but which describe only a local 
region of the system of interest. The proper concept that 
encodes the similarity between such mixed states is that 
of the reduced fidelity, which is the maximum pure state 
overlap between purifications of the mixed states Q. ft 
has proven useful in the analysis of a number of ordi- 
nary symmetry- breaking quantum phase transitions Q. 
But since the reduced fidelity is a local property of the 
system, similarly to that of a local order parameter, one 
may think that it would be less sensitive to a TQPT, 
which involves a global rearrangement of nonlocal quan- 
tum correlations However, this intuition turns out to 
be wrong. As we show in this paper, several TQPTs are 
accurately signaled by a singularity in the second-order 
derivative of the reduced fidelity. Moreover, the singular- 



ity can be even stronger than for the (pure state) global 
fidelity. The fact that a TQPT gets imprinted in a local 
quantity may at first seem surprising, but, as we shall see, 
parallels and extends results from earlier studies [l^, [ll[ . 

Fidelity and fidelity susceptibility — The fidelity 
F[(3, f3') between two states described by the density ma- 
trices p{/3) and p{l3') is defined as [1] 



F{f3, f3') = Tl■^/7Mm)7M■ (1) 

When a system is in a pure state, = 
F{P,f3') becomes just the state overlap | (*(/3')|*(/3)) |. 
When the states under consideration describe a subsys- 
tem, they will generally be mixed states, and we call the 
fidelity between such states reduced fidelity. In the limit 
where (3 and f3' = f3 + 5(3 are very close, it is useful to 
define the fidelity susceptibility jl2| 
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consistent with the pure state expansion Fk, I—xf^P"^/"^- 
The Castelnovo-Chamon Model ~ The first model we 



consider was introduced by Castelnovo and Chamon IC | , 
and is a deformation of the Kitaev toric code model 13 1 . 
The Hamiltonian for N spin- 1/2 particles on the bonds 
of a square lattice with periodic boundary conditions is 



if = -Ao ^ Bp - Ai ^ A, + Ai ^ e-'3 5:>. 



(3) 



where As = HiGs ^"^"-^ ~ IliGp ^1 ^^^'^ 
and plaquette operators of the original Kitaev toric code 
model. The star operator As acts on the spins around 
the vertex s, and the plaquette operator acts on the 
spins on the boundary of the plaquette p. For Ao,i > 
the ground state in the topological sector containing the 
fully magnetized state |0) is given by [l3| 
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with Z{I3) = Egece'^^''''^^^ where G is the Abelian 
group generated by the star operators As, and al{g) is 
the z component of the spin at site i in the state 5|0). 
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FIG. 1; (Color online.) Mapping between the Castelnovo- 
Chamon model and the 2D Ising model. The spins of the 
former reside on the lattice bonds (filled black circles), and 
the spins of the latter on the vertices. Left: af = OsOgi , where 
i is the bond between the neighboring vertices {s,s'). Mid- 
dle and right: For i and j nearest (next-nearest) neighbors, 
the mapping gives {afaj) — {OsOgiOgiiOs) = {ds'Ss"), where 
(s',s") are next-nearest (third-nearest) neighbors. 



When /? = the state in ([4]) reduces to the topologically 
ordered ground state of the toric code model . When 
/9 — > oo the ground state ^ becomes the magnetically 
ordered state |0). At Pc = (l/2)ln(V2 4- 1) there is a 
second-order TQPT where the topological entanglement 



entropy Stono goes from Stopo = 1 for (3 < f3c to Si 
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for P > Pc [10|. The global fidelity susceptibility xf close 
to Pc was obtained in Ref. 0, and found to diverge as 



XF-ln\Pc/P-l\. 



(5) 



We here calculate the single-site reduced fidelity be- 
tween the ground states of a single spin at two different 
parameter values P and P' . To construct the density ma- 
trix Pi for the spin at site i we use the expansion 



(6) 



/j=0 



with = Ij, and with the expectation values taken 
with respect to the ground state in There is a 

one-to-two mapping between the configurations {g} = 
G and the configurations {9} = of the classical 
2D Ising model H = -J E<s.s'> ^s^'s' with 6^ ^ 
— 1 (+1) when the corresponding star operator As is 
(is not) acting on the site s Thus ct| = 9s9si, 

where i is the bond between the neighboring vertices 
{s,s'), see Fig. [H This gives {GSiP)\af\GS{P)) = 
{l/Z{P))j:,^^9s9s'e^^<^"y")'^"'^"' = E{P)/N, where 
P is identified as the reduced nearest-neighbor cou- 
pling J/T = P of the Ising model with energy E{P). 
The two expectation values (G5'(/3)|d-f |G'S'(/3)) and 
{GS{P)\a\\GS{P)) are both zero, since {Q\gdfg'\Q) = 0, 
V(7,(7' e G, and similarly for . It follows that pi = 
{l/2)<Xia.g{\ + E{P)/N, I- E{P)/N) in the af eigenba- 
sis. Since the density matrices at different parameter 
values P and P' commute, the reduced fidelity ([1]) is 
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FIG. 2: (Color online.) Single-site fidelity (a), single-site 
fidelity susceptibility (b), two-site fidelity (c) and two-site 
fidelity susceptibility (d) of the Castelnovo-Chamon model 
calculated with a parameter difference 5/3 = 0.001 and with 
N ^ oo. The reduced fidelity susceptibilities will diverge ac- 
cording to Eq. (|9]) when 5(5 — » 0. In (c) and (d) we plot for 
both nearest (NN) and next-nearest (NNN) neighbors. 



where {A^} ({A^}) are the eigenvalues of pi{P) {pi{P')). 
The energy E{P) of the 2D Ising model in the ther- 
modynamic limit iV — > oo is given by E{P)/N = 
-coth(2/3) [1 + (2/7r)(2tanh2(2/?) - 1)K{k)\ /2, where 
K(n) = Jo^^d9{l - K2sin2 6»)-i/2 and k = 
2sinh(2^) /cosh^(2/3) This gives us the plot of the 
single-site fidelity shown in Fig. [51 where we see that the 
TQPT at Pc =^ (l/2)ln(V2 1) « 0.44 is marked by a 
sudden drop in the fidelity. 

The single-site fidelity susceptibility xf is 
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for commuting density matrices |la |. Here 813X12 = 
±{2N)-'^di3E{P)=±{2NP^)-^C{P), with C{P) the spe- 
cific heat of the 2D Ising model. Thus xf diverges as 



XF-ln\Pc/P-ll 



(9) 



at Pc, faster than the global fidelity susceptibility in ([5|). 
In Fig. [2] we plot the single-site fidelity susceptibility 
using Eq. but with finite SP = 0.001. 

The two-site fidelity can be obtained in a similar way. 
We expand the reduced density matrix pij as 



(10) 



The only non-zero expectation values in (fTO|) are 
(aOaO) = 1, (afaO) = (af), (aOa|) ^ and (afap. 
Translational invariance implies that (ct^) = (irf), so that 



k,=-il + {at){at+a^) + {ala^)ala^). 



(11) 



Fip, P') = i:rJpAp)pdp') = V VA^, (7) The eigenvalues are seen to be Ai,2 = (1/4)(1 ± 2((Tf ) + 

i (af<7|)) and A3,4 = (1/4)(1 - (<Tfa|)). Now the fidelity 
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can be calculated using Eq. Here we focus on the 

cases where i and j are nearest and next-nearest neigh- 
bors. Then the mapping to the 2D Ising model gives that 
(af) = {0s9s') and (o-J) = {9,," 9s'"), where i (j) is the 
bond between the neighboring vertices (s, s') ((s",s"')). 
When i and j are nearest (next-nearest) neighbors, we 
get (afa|) {9s9.,'9s"9s) = {9s'9s"), where {s',s") are 
next-nearest (third-nearest) neighbors on the square lat- 
tice (cf. Fig.[l|). As before, (af) = E{P)/N. We obtain 
{9s' 9s") from the equivalence between the 2D Ising model 
and the quantum ID XY model 
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+ M^), (12) 



where a± = (1 ± 7)/2. This has been shown to give 

{9oS)9n,n) = (o'o'3'5^)xy|7 = l,/t=(sinh2/3)-2 

for Ising spins on the same diagonal, and 

" ^ -^cosh^{(i*){al^ai^,)xY\^=^p,h=hp 
-sinh\n{^l,^l'): 



(13) 
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for Ising spins on the same row (or, by symmetry, col- 
umn), where tanh/3* = e~^'^, 7^ = (cosh 2/3*)"-^ and 
hfj = (1 - 72)i/2/tanh2/3 [11]. Known results for the ID 
XY model give [ill 
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where G,-' = (I/tt) d(j){h — cos(j))cos{(j)r')/A^{h) + 
(j/tt) d(f) sin (f)ain{(l)r')/Aci,{h) and A^{h) ~ 
((7sin0)^ + (/i — cos t/))^)"'^/^. These relations allow 
us to plot the two-site fidelity, and also the two-site 
fidelity susceptibility using Eq. ([2|), see Fig. [2] Note that 
the two-site functions are only slightly different depend- 
ing on whether the two sites are nearest neighbors or 
next-nearest neighbors. It follows from Eq. ([8]) that also 
the two-site xf has a stronger divergence at criticality 
than the global fidelity susceptibility. 

It is interesting to note the slight asymmetry of the 
reduced fidelities around the critical point, seen in Fig.[2l 
indicating a somewhat smaller response to changes in 
the driving parameter in the topological phase. 

The transverse Wen-plaquette model — We now turn to 
the transverse Wen-plaquette model, obtained from the 
ordinary Wen-plaquette model [l3| for spin- 1/2 particles 
on the vertices of a square lattice by adding a magnetic 
field h [H, 



(17) 



where Fi = af a^^-a^_^_^_^_-a^_^_. and g <0. The boundary 
conditions are periodic. At /i = the ground state is the 
topologically ordered ground state of the Wen-plaquette 
model jl8| and in the limit — > 00 the ground state is 
magnetically ordered. Since Fi, ctJ have the same com- 
mutation relations as T[_^^^2+y/2' ^j-i/2+y/2^f+i/2-y/2 

(where the Pauli matrices f act on spin- 1/2 particles at 
the centers of the plaquettes), the Hamiltonian (|17p can 
be mapped onto independent quantum Ising chains. 



H = -h 
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(18) 



with gj = g/h^ and where r^.i and are the im- 

ages of afa^_^_~af^^J^~a^_^_~ and af respectively ^M- The 
index a denotes the diagonal chains over the plaquette- 
ccntered sites, and i is the site index on each diagonal 
chain. Known results for criticality in the quantum Ising 
chain imply that the transverse Wen-plaquette model has 
a TQPT at g/h = 1 [H. 

We now calculate the reduced fidelity. The mapping 
onto the quantum Ising chains immediately gives that 
i^i) = (T^i'^^^i)- the af basis, the Hamiltonian 
(fTT]) only flips spins in pairs, therefore we get (af) — 
and (af) = 0. The single-site reduced density matrix 

© is therefore given by = (1/2)(1 + 1 Tf, 1 )CTf ), 

2 ^ 2 

which is diagonal in the af basis, with eigenvalues A1.2 = 

(1/2)(1 ± (T^_i,Tf^)). The single-site fidelity is thus 

* 2 2 

given by Eq. Q, and (f^_ 1 1 ) is calculated using Eq. 

^ 2 2 

p5|) with 7 = 1 and h ^ gj. The result reveals that the 
TQPT is accompanied by a sudden drop in the single- 
site fidelity. Now. dg,\i^2 = '^\dg,{T^_xT^^i) , which di- 
verges logarithmically at the critical point = 1. There- 
fore Eq. ([5]) implies that at h/ g = 1, xf diverges as 



XF^ln' \g/h-l\. 



(19) 



as in Eq. ([9]) for the Castelnovo-Chamon model. 

We can also calculate the two-site fidelity for two near- 
est neighbor spins at sites i,]. All non-trivial expecta- 
tion values in the expansion ()10p of the reduced den- 
sity matrix, except {af), (ctJ) and {af&j), will be zero, 
since only these operators can be constructed from those 
in the Hamiltonian (fT7|) . The mapping onto the quan- 
tum Ising chains gives {afaf) = {T^_iTf,iT^_if^',i) = 

J * 2 2 ^ 2 ' 2 

((ri*^ j_T?^ 1 ))^. Thus the two-site density matrix is 
given by' p,, ^ (1/4)(1 + (ff_ ^ ff^ ^ )(af + aj) + 
{{^'i-.i^'i+i))'^^i'^^); which is diagonal in the afafj eigen- 
basis. The eigenvalues arc Ai,2 = (1/4)(1 ± {t^_iT^_^i))^ 
and A3.4 = (1/4)(1 - ((Ti"_if:"_^i))2). Taking deriva- 
fives of the eigenvalues Ai^2,3,4 and inserting them into 
Eq. ([8]) shows that also the two-site fidelity susceptibility 
diverges as xf ^ In^ \g/h—l \ ath/g = 1. Contrary to the 
case of the Castelnovo-Chamon model, xf for one and 
two spins now diverges slower than the global fidelity sus- 



4 



ccptibility, which shows the Yp ^ \g/h—\\~^ divergence 
of the quantum Ising chain [20|. 

The Kitaev toric code model in a magnetic field — 
Addin g a magnetic field h to the Kitaev toric code 
model [Tsj gives the Hamiltonian [ll| 

H = -XoY,Bp-XiJ2^s-hY,'^t, (20) 

p si 

where the operators Bp and Ag are the same as in Eq. ([3]). 
In the limit Ai ^ Aq, h, the ground state \GS) will obey 
As \GS) = \GS), Vs. Then there is a mapping to spin- 1/2 
operators f acting on spins at the centers of the plaquc- 
ttes, according to Bp i-^ t^, af ^ Tpf^. Here i is the site 
shared by the two adjacent plaquettes {p,q). This maps 
the Hamiltonian ([201) onto [Hi 

i^ = -Ao^f;-/^5]f^%^ (21) 
p (p-q) 

which is the 2D transverse field Ising model with mag- 
netic field Xo/h = h' . Now. the mapping tells us that 
(af) = {TpT^), and the symmetries of the Hamilto- 
nian ^ imply (o-f) = and (af) = 0. The single- 
site reduced density matrix is therefore given by pi = 
(1/2)(1 -I- {Tpfg)af), which has the same form as in the 
transverse Wen-plaquette model. Since numerical results 
have shown a kink in (Tpf^) at the phase transition at 
/i^ « 3 [ll|, it follows that the single-site fidelity will 
have a drop at this point. Further, the divergence of 
dh'{TpTg) at the critical point implies a divergence of the 
single-site fidelity susceptibility at h'^. Thus, the scenario 
that emerges is similar to those for the models above. 

Discussion — To summarize, we have analyzed the re- 
duced fidelity at several lattice system TQPTs and found 



that it serves as an accurate marker of the transitions. In 
the case of the Castelnovo-Chamon model [l^] , the diver- 
gence of the reduced fidelity susceptibility at criticality 
can explicitly be shown to be even stronger than that 
of the global fidelity @. Our analytical results rely on 
exact mappings of the TQPTs onto ordinary symmetry- 
breaking phase transitions. Other lattice models exhibit- 
ing TQPTs have also been shown to be dual to spin- 
[2]| or vertex [H] models with classical order, suggesting 
that our line of approach may be applicable also in these 
cases, and that the property that a reduced fidelity can 
detect a TQPT may in fact be generic. While counter- 
intuitive, considering that the reduced fidelity is a local 
probe of the topologically ordered phase, related results 
have been reported in previous studies. Specifically, in 
Refs. [l^ and the authors found that the local mag- 
netization in the Castelnovo-Chamon model and the Ki- 
taev toric code model in a magnetic field, while being 
continuous and non-vanishing across the transition out 
of topological order, has a singularity in its first deriva- 
tive. The fact that local quantities can spot a TQPT 
is conceptually satisfying, as any physical observable is 
local in nature. Interesting open questions are here how 
the concept of reduced fidelity can be applied to TQPTs 
in more realistic systems, such as the fractional quan- 
tum Hall liquids, and how reduced fidelity susceptibility 
singularities depend on different topological and classical 
orders involved in the transitions. 
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